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ABSTRACT 



The subject of this paper is an investigation of the 
significant factors that affect the performance of a sub- 
cavitating pitot type intake for a waterjet propulsion 
system for hydrofoil craft. A computer model for design 
of an inlet nacelle is developed which accepts values of 
two parameters and several system constants from a control 
program and designs a nacelle to meet those requirements 
subject to limitations imposed by cavitation and geometry. 

The parameters are inlet velocity ratio and inlet diameter 
to maximum diameter ratio. A body of revolution config- 
uration with profile shape based on certain desirable 
pressure distributions is assumed. Extensive use is made 
of empirical data from the available unclassified literature. 

The design conditions are those existing at maximum 
speed, called cruise speed, and the minimum foil-borne 
speed, called take-off speed. 

The subroutine will accept up to three additional speeds 
at which a performance prediction will be made after design 
is completed from th first two. 

Pressure recovery characteristics and cavitation limits 
for some typical designs are displayed as graphs. 

The subroutine listing and flow diagram is included in 
Appendix C. In addition, an alternate design method which 
accepts only the inlet velocity ratio as a parameter is 
included in Appendix D. 
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Chapter 1. 

PROBLEM STATEMENT 

Ram (pitot) type inlets for water jet propulsion systems 
in hydrofoil craft are responsible for a significant part 
of both the external drag and internal flow losses and 
have a strong influence on the propulsive efficiency of 
the system. 

The objective of this investigation is to identify and 
quantify the significant characteristics of a ’’well de- 
signed” inlet in a manner that will allow their use at the 
preliminary design stage to assist in rapidly determining 
the feasibility and relative merit of various propulsion 
system configurations in a craft of a given size with 
given speed and range requirements. 

The term, ’’well designed”, as used with reference to 
the inlet and nacelle will mean that inlet configuration 
which best suits the propulsion system as a whole when 
optimized with respect to some chosen criterion. 

A computer subroutine is to be developed for use with 
an optimization program concurrently being developed by 
R. P. Gill (ref. l) which calls this subroutine and others 
written by Gill and R. C. Percival (ref. 2). The optimizat- 
ion criterion chosen for this particular program is least 
propulsion system weight including fuel, but the inlet 
subroutine, called SUBROUTINE NAGEL, does not depend on 
that particular criterion and could be used for others. 

Engineering approximations and empirical fits are 
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used wherever they are suitable and consistant with the 
purpose of the study. 

There are at least three broad categories of nacelles 
which could be used for ram inlets; subcavitating, super- 
cavitating, and base-vented. Only subcavitating nacelles 
are considered in this study. Persons who are interested 
in supercavitating and base-vented nacelles may wish to 
read references (3) and (4) for further discussion on 
those types. 

Subsonic aerodynamics provides much of the data 
utilized in this study, but hydrodynamic design has the 
additional complication of cavitation avoidance. For 
the purpose of this paper cavitation is assumed to start 
when the local static pressure decreases to the vapor 
pressure of sea water. Incipient cavitation is imposed 
as an absolute limit on the decrease in static pressure 
at the design conditions even though a certain amount of 
cavitation may be allowable in a real life situation. 
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Chapter 2. 

DESIGN CONDITIONS 

There are three possible flow conditions at the inlet 
and Bernoulli's equation indicates the changes in local 
static pressure coefficient associated with each of them. 

When Vi/Vo *^1, the stream lines are diverging (Fig.l) 
and since there are no losses (to the inlet plane, at 
least), local velocities inside are lower than free stream 
velocity causing the pressure coefficient to be positive 
and no cavitation danger exists. Local velocities on the 
exterior of the forebody are higher than in the free stream 
and the pressure coefficients are negative (Fig. 2), so 
cavitation danger exists and increases as Vi/Vo decreases. 
Thus, for any particular nacelle and free stream velocity 
combination, there is a minimum velocity ratio below which 
cavitation occurs in the external flow. 

When Vi/Vo-^1, the stream lines are converging 

(Fig. 1). Local velocities in the external flow will 

still be somewhat greater than in the free stream but 

lower than in the diverging flow case, so if the external 

flow was cavitation free when Vi/Vo ^1, it will also be 

cavitation free when Vi/Vo^l at the same Vo. Local 

velocities in the internal flow are now greater than the 

free stream with a corresponding decrease in static pressure 

coefficient, and cavitation inside the lip is a definite 

possibility. There is therefore a maximum velocity ratio 

for each nacelle and Vo above which cavitation will occur 
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Figure 1. Entering stream tube illustrating two flow 
conditions into the nacelle. 
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Figure 2 * Minimum Cp in the external flow vs . for 

several of the Lockheed nacelles 
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in the internal flow. 

When Vi/Vo = 1, the flow is in an intermediate state 
between the other two conditions and no new limitations 
are imposed. This condition is of no interest, then, from 
the design viewpoint. 

Changes in local velocity due to the geometric con- 
figuration of the nacelle, lip shape, bluntness of the 
forebody, etc., are perturbations on the flow conditions 
described above and are assumed to be minimized by the 
detailed design of the nacelle. No attempt is made here 
to define that detailed design except to assume that it is 
at least as good as the best of the empirical examples 
and data available to this investigation. 

A very frequently encountered set of requirements, 
as indicated by the literature, includes the following 
situation: 

Cruise speed is approximately twice take-off speed. 

Cruise flow rate is approximately equal to take-off 

flow rate. 

Let the subscripts 1 and 2 designate cruise and take-off 
conditions, respectively. If the inlet area is a constant, 

Ql=Vi]L Ai 

Vii=Vi2 

but Voi=Vo2*2, so Vi/Voi=]^Vi/Vo 2 . 

If Vi/Voj^^O.7 to 0.^ as is frequently the case, 

(refs. 4,5,6,) then Vi/Vo2=1.4 to 1.6, but except when 

the inlet is very deeply submerged, these high velocity 
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ratios will cause cavitation inside the inlet, incurring 
high losses and perhaps choking the inlet. 

One way of decreasing this effect is to decrease the 
inlet velocity ratio at cruise by increasing the area. 

This approach very quickly reaches the minimum velocity 
ratio below which cavitation occurs in the external flow 
with its accompanying high drag and possible material 
damage . 

As speed increases, the interval between these two 
limits becomes smaller, so the problem grows increasingly 
difficult as attempts are made to achieve higher speeds. 

A solution to the problem may be found in the use of 
auxiliary inlets to increase the effective inlet area when 
high velocity ratios are required. 

If the existance of an auxiliary inlet is postulated, 
the velocity ratio at take-off is no longer coupled to 
the velocity ratio at cruise and both become design con- 
ditions as the total inlet area for each must be determined 
independently. 

An additional limitation is imposed by the elbow 

where the internal flow leaves the nacelle and enters the 

strut. This elbow will normally contain turning vanes to 

help minimize losses. Since the local velocities on the 

turning vane surface will be greater than the average 

velocity in the flow entering the elbow, a danger of 

cavitation exists. This turning vane cavitation is avoided 

by diffusing the flow in the nacelle to reduce the 
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velocities and increase the static pressure. Therefore, 
a minimum diffusion ratio dictated by a maximum nacelle 
exit velocity, above which cavitation occurs in the elbow, 
is a third design limit. This limit must be checked at 
both cruise and take-off because it is a function of both 
the flow rate and the losses. 

In summary, then, both cruise and take-off are design 
conditions and cavitation avoidance establishes absolute 
limits on the minimum velocity ratio at cruise, the maximum 
velocity ratios at both cruise and take-off, and the 
minimum diffusion ratio. 
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Chapter 3» 
NACELLE SHAPE 



The principles of design used in jet engine nacelles 
and air intakes, already established in the field of 
aerodynamics, can be used for guidance in the selection 
of water jet inlet nacelles (ref. 3 and 4) • 

Kiichemann and Weber (ref. 7) have developed the de- 
finition of a theoretical optimum nacelle contour by con- 
sidering the momentum balance between the free stream and 
the entry plane of a semi-infinite intake. The resultant 
force on the nacelle is given by 

V i '"i 

but since there can be no resultant force on a semi-infinite 
body in potential flow, this force must be equal and 
opposite to the integral of the pressure over the forebody 
area of the nacelle. 



^ ^m 

V J % dA = -4/’v2 

h Ai 

If the pressure coefficient is a constant over the 
entire forebody then the local velocity, V^, is a constant 
and the projected frontal area, required to balance 

Fjj is defined by the above equation. When the local 
velocity is limited to some maximum value to satisfy Mach 
number requirements or, in the case of water jet nacelles, 
cavitation requirements, then the minimum projected frontal 



m 



[l - (Vj^/V^)2] dA 
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area and will satisfy the relationship 



Am & 1 + (1 - Vi/Vo)2 

Ai (Vmax/Vo)2 - 1 



The optimum intake contour of Kiichemann and Weber is 
one which at its design achieves the goal of keeping 

Vmax constant and at its specified value over the entire 
forebody and the equals sign applies. 

A^/A^ can be represented by the square root of its 
inverse, in the case of circular intakes. 

(1 - Vi/V„)2 



ri/Dm - 



1 + 



(V /v«)^-l 
' max' 

It can be seen that a maximum value of D./D_ is defined 

1' m 

by a minimum value of for any given V^. One of the 

design limits chosen in the previous chapter was avoidance 
of cavitation in the external flow at the relatively low 
velocity ratio at cruise speed. This can now be translated 
into a maximum value of D^/Djjj. 

It is not practical to build such perfect intakes 
although intakes with reasonably uniform pressure dis- 
tributions at certain velocity ratios have been developed 
(ref. 3, 4, 7). In general, the velocity distribution 
will have peaks. 

The generation of practical contours with acceptiable 
pressure distributions has been largly a matter of 
experimental trial and error. Kiichemann and Weber de- 
scribe three sets of geometrically similar shapes which 
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they designate as Classes A, B, and C» The geometric 
description of these are given in Table 1 and Figure 3 • 

References 3 and 4 describe two sets of intakes 
whose pressure distributions were generated by numerical 
solutions to the potential flow problem. The Lockheed in- 
takes were chosen for input data for this study chiefly 
because the data covered a wide range of geometric variation 
and data on the interal flow distribution was available 
for each of the individual inlets. The major geometric 
parameters and peak pressure coefficients for each of the 
nine inlets is given in Tables 2 through 11. The external 
contour of these nacelles is very close to, but not quite, 
elliptical. The internal lip contour is elliptical. The 
Kiichemann and Weber intakes are elliptical both internally 
and externally. Brown and Traksel state that this 
variation has an influence that is secondary to the length 
to diameter ratio, L/D^^^, and Figure 4 shows the 

similarity of the major parameters. The dashed line was 
tabulated as the basic geometric description of nacelle 
forebodies for SUBROUTINE NAGEL. The three data points 
to the left represent intakes which had lower pressure 
coefficients than the others, so the dashed line represents 
the best of the available data points . 

To obtain a functional relationship for the variation 
of pressure coefficient with geometery, the peak pressure 
coefficient at each velocity ratio and at each angle of 
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attack was plotted against the value of and also 

against L/Djj^. L/Djj^ was chosen as the appropriate parameter 
to use because it gave a smoother distribution of data 
points. The desired relationship was then obtained by 
drawing a smooth envelope that enclosed all points and 
passed through the points of highest C^. This envelope 
is considered to represent the "well designed" inlet 
mentioned in the problem statement based on the assumption 
that if a nacelle has already been designed which falls 
on that line (the discrete data points through which the 
line was drawn), then a nacelle could be designed for any 
other point on that line. The plots are shown in Figures 5 
through 10. Figures 11 through 16 show similar plots for 
the internal flow conditions. 

There can be no losses in the stream tube from out 

in the free stream up to the inlet plane, but from the inlet 

plane on into the ducting of the nacelle friction and 

turbulence take a toll of the head available in the internal 

flow. The amount of head loss is primarily a function of 

velocity ratio, Reynolds number and lip shape. Blackaby 

and V/atson investigated experimentally several lip shapes 

including circular profiles, elliptical profiles and sharp 

profiles. In general, the elliptical profile caused less 

losses at a given Vj^/V^ than the others. The particular 

results for the lip which they designated ISE is given as 

a function of Reynolds number for several velocity ratios 

in reference 3, Figure 74 of Report No. 5. The loss 
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coefficient, defined as AP/q^, approaches a constant value 
for each velocity ratio for Reynolds numbers greater than 
^xlO^. These constant values as read from the citation 
above are given in Table 12. The Ram Pressure Recovery, 
defined as 

RPR= (P - Po)/(Po - Po) = 1 - AP/q„ 
is plotted in Figure 17, and points from that plot are 
tabulated in the subroutine. 
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Chapter 4 
THE DIFFUSER 

The inlet nacelle modeled in this study is assumed to 
consist of a streamlined body of revolution containing 
a nose intake opening, a diffuser and an auxiliary inlet. 
The forebody and intake lip serve to guide the flow into 
the internal ducting. That internal ducting includes the 
diffuser which has at least two primary functions. One 
is to slow the flow velocity to below the critical value 
for cavitation in the strut elbow. The other is to de- 
crease the flow velocity further to reduce ducting losses 
if it benefits the system as a whole. 

Three main profile shapes of diffusers with circular 
cross section can be considered. They are bell, trumpet 
and straight walled conical diffusers. To avoid confusion 
of terms, a bell diffuser is one having its mean diameter 
greater than the mean diameter of a straight walled conical 
diffuser with the same inlet and exit areas while a trumpet 
diffuser has a mean diameter less than that of the cone. 

At similar inlet flow conditions, the bell shape results 
in a lower total head loss than the other two and the 
trumpet results in a higher loss {ref. 4, 9)» 

In general, the diffuser with the least loss in total 
head is preferred for this application if other con- 
siderations are equal because the difference in losses must 
eventually be made up by expenditure of prime mover power. 
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Current knov/ledge of boundary layer growth and be- 
havior in adverse pressure gradients does not allow accur- 
ate prediction of diffuser performance by analytic methods 
(ref. 3), and very little empirical data on bell diffuser 
performance that is directly applicable was found. Ref- 
erence 9 contains an experimental comparison of the three 
shapes which indicates that conical diffusers have pressure 
recovery characteristics very close to the bell and sig- 
nificantly better than the trumpet. Experimental results 
contained in references 9 through 13 can be applied to 
predict the losses of straight walled conical diffusers 
and can be considered a lower bound of that which can be 
expected from a bell. 

The total pressure loss for a conical diffuser can 
be predicted by 

AP = (C 4 f L/djj^) qi, (ref. 13), 
where the second term is a ducting loss of a duct the same 
length as the diffuser and of a diameter equal to the mean 
diameter of the diffuser. For a cone, 

d„ = (D^ 4 D2>/2 . 

The first term accounts for diffusion losses as a percent- 
age of the losses in a sudden expansion. is the sudden 

expansion loss factor and is given by 

= (1 - Ai/A2)^ . 

The proportionality factor C is a function of ex- 
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i 



1 



pansion angle, 9, or the double angle, 20, and is determined 
experimentally as in reference 13 • 

The minimum loss for any given diffusion ratio occurs 
when the value of 20 is in the vicinity of 5° to 10°. Flow 
separation starts interfering with the flow and partially 
stalling the exit when 20 is in the vicinity of 20° 

(ref. ll), therefore the region of interest is for 20 less 
than 20°. In this region the curve of C vs. 20 can be 
approximated by 

C = 3.190 X 10-3 (0)2 + x 10“^ (0) 

where 0 is in degrees. 

There appears to be no reason in this application to 
decrease the angle below the value of minimum loss because 
the diffuser length increases and friction losses start to 
rise accordingly. For the computer model limits of 3° and 
10° were set on the minimum and maximum values of 0 which 
has the effect of establishing a maximum and minimum length. 

The desired total length of the nacelle is generally 
established by external hydrodynamic flow, structural and 
other considerations that do not include the diffuser and 
other internal length requirements. If the internal com- 
ponent lengths are less than the desired length, they will 
all fit within the envelope. If not, additional external 
drag due to the added length and wetted surface is generated 
and the power required to overcome that drag increment can 

be charged to the diffuser length since other internal 
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lengths are fixed. This becomes a criterion for determin- 
ing the .proper diffuser length. The sum of the power re- 
quired to overcome diffuser losses and the additional 
power due to the increment in drag is minimized with re- 
spect to diffuser length. The calculations used in this 
determination are given in Appendix B. 

A limit on the minimum diffusion ratio is given by the 
requirement that the maximum diffuser exit velocity must 
be less than the critical velocity which will cause 
cavitation in the elbow. Experimental results in reference 
11 give exit velocity profiles for six diffuser examples, 
in two of which 20 = 10® and in four of which 20 = 12°. 

The diffuser area ratio was 2.0 in all cases, entry Mach 
number was 0.26, and boundary layer thickness ranged from 
0.003 to 0.03 times the entry radius. Integration of 
these distributions graphically gives a range of 
^max./^ave. ” to 1.50, with a mean value of 1.33 

over the six cases. However, since boundary layer thick- 
ness has poor predictability, the average exit velocity is 
not known but it will be somewhat greater than the flow 
rate divided by the exit area. If V is calculated by 

^max. ^ ^ Q/^2» 

k = 1.33 will probably give a low value and k = 1.50 may 
give a high value but will be on the safe side. 

The upper limit on will be the critical value 

of V 2 determined from Appendix A. 
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Chapter 5» 

THE AUXILIARY INLET 

It was indicated in chapter 2 that an auxiliary inlet 
was probably required to obtain the flow rates necessary 
at take-off speed and remain subcavitating at cruise. 

The possible configurations for providing auxiliary 
inlets are many and varied. Slots in the nacelle wall with 
doors that open inward, leading edge openings on the strut 
and a movable nacelle lip are all possibilities, to mention 
only a few. 

The model for this study was assumed to be of the first 
type. Figure 3 shows the axixiliary inlet as dashed lines 
in the profile. It is assumed that the inlet is in the 
form of a series of openings in an annular ring around 
the nose opening with some fraction of that ring being 
structure to support the nose. The dimension X in the 
figure is^ the width of that annular ring in the direction 
normal to the inlet flow through the opening. 

The auxiliary inlet flow is assumed to enter and 
combine with the flow from the nose opening prior to enter- 
ing the diffuser because that is the area of lowest static 
pressure and, since the inlet will only be used at take-off 
when the flow rate is high, it is important that the entire 
combined flow pass through the diffuser to prevent 
cavitation in the elbow. 

When the inlet is closed the length of ducting, L , 
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associated v/ith it is considered to be a short length of 
pipe and the loss for that length of pipe is added to the 
other system losses. 

This type of inlet appears to be in much the same flow 
conditions as the ’’gliding scoop” described in reference 14« 
Since the auxiliary inlet is open only when the system de- 
mands a high flow rate its inlet velocity ratio is estimat- 
ed to be about 1.0 or slightly higher. In this condition 
the lip loss coefficient is approximately 0.2 (ref. 14) • 

This value is the one used in the subroutine for auxiliary 
inlet losses. 



- 29 - 



Chapter 6. 

DESIGN RATIONALE 

The design of an inlet for the specific conditions 
imposed by the optimization program is based on comparison 
of the characteristics implied by the system parameter 
values with those stored in the data tables and sizing the 
individual components to remain cavitation free at the 
design conditions. The subroutine listing and a computer 
generated flow diagram are in Appendix C. This and the 
following chapter will attempt to clarify the reasoning 
used in its derivation. 

When the program calls SUBROUTINE NACEL, all physical 
constants, craft configuration, depth of submergence, 
free stream velocity, flow rate and angle of attack for 
each design point are known. Inlet velocity at cruise 
and are passed as system parameters. Actually, 

''lAo at cruise is the system parameter, but it is cal- 
culated from and rather than passed. 

Initially, free stream static pressure, free stream 
dynamic pressure, incipient cavitation number, and total 
pressure for both cruise and take-off speeds are calculated. 

<3o = 

<Toi = (Po - Pv>i/’'^o 

Po = Po * % 

The length to diameter ratio, L/D^, of the nacelle for 
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v/hich the external minimum pressure coefficient is equal 
to the negative of the free stream incipient cavitation 
number at cruise is determined by interpolation in the 
data table where 



L/D = f(C , V./V , Cl). 

' m ' p» i' o* 

The inlet to maximum diameter ratio, correspond- 

ing to that L/Dj^ is found in the tabulated geometric 
characteristics. The value obtained defines a nacelle for 
which the external flow on the surface is at the point of 
incipient cavitation and therefore is a maximum value. 



(D,./D„) = f(L/D„) 

If the Di/Djjj passed by the calling program is greater 
than (^i/^m^max* particular set of parameters is 

rejected by setting the nacelle weight to a very high 
number and returning to the calling program. If 0^/0^^ is 
less than then it is cavitation free and the 

design proceeds. 

Inlet area is determined from the flow rate and inlet 
velocity at cruise which then allows D. and D to be 
determined. 

The maximum allowed velocity ratio for both speeds is 
determined by interpolation in the data table again. 

(V./V ) = f(L/D , cr . , ot) 

This limit is imposed by cavitation -on the inside of the lip. 

If at cruise, V^/V. is greater than (V./V ) , the 

parameters are rejected in the same manner as before. If 
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not, the take-off condition is checked to determine if 
an auxiliary inlet is required. The maximum flow rate 
allowed by (V./V^) enters through the nose inlet, the 
remainder must enter through the auxiliary inlet. 

The pressure recovery of the lip is a tabulated 
function of so the total pressure inside the inlet is 

h ‘=0 Po 

and the static pressure inside is given by 



Pi ^ Pi - . 

Pressure recovery of the auxiliary inlet is known 
from chapter 5 so the total pressure inside the auxiliary 
inlet is 



^aux “ ^^^aux *^o ‘‘‘ Po * 

The static pressure at a point is unique so the static 
pressure inside the auxiliary inlet is the same as that 
previously calculated from the flow inside the nose inlet. 
Therefore, the velocity through the auxiliary inlet is 

''aux ■- ^ • 

The area of the auxiliary inlet is found by 

^aux " "^aux/ ^aux • 

The total pressure of the. combined flow is estimated 
as the mass weighted average of the combining flows. 

^c “ (*^aux Paux ■•■ •Qi PfV Qq 
where ^ ^i • 

The average velocity of the combined flow is 
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calculated from the dynamic pressure 

Vi = [ 2(p^ - p.)//>] i . 

and the pressure recovery coefficient to this point is 

RPRc = <'’c - Po>/ % • 

The diffuser entry and exit diameters are set by 

^1 ■ °i » °2 ^ Djjj . 

The 0.9 factor is somewhat arbitrary. Some fraction 
of the nacelle diameter must be allowed for structural 
requirements and it was felt that ten per cent would be 
sufficient in most cases. 

The internal component length requirements are 
calculated in the manner of Appendix B and the sum of 
those lengths plus the lip length is used as the nacelle 
length subject to a minimum value of 5»5 0^^^* This factor 
is also arbitrary and is the average of the recommended 
values of found in the literature. Hoerner (ref. 15) 

shows from empirical data that the minimum drag nacelle 
on a wing has a value of approximately 3»0, but since 

the internal length requirements are always greater than 
3.0 (see Appendix B), it was felt that value might place 
an unncesssarily high drag penalty on the diffuser. 

The larger of the two flow rates, cruise or take-off, 
is used for sizing the diffuser. 

The diffuser loss, lip loss, and pipe loss, if any, 

are calculated for both cruise and take-off. Their sum in 

each case is used in the equation derived in Appendix A 
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to determine the critical velocity value at the diffuser 
exit for cavitation on the elbow turning vane. 



'^crlt = [‘^ol ^ 1 - Ploss/i/’''o ] ^ ^ ^tvi)^ 

The maximum diffuser exit velocity as determined in 
chapter 4 is 

''a.ax = "-5 • 

These two velocities are compared for both conditions. 

If V is greater than V for either condition, the 

crxo 

subroutine returns to the calling program as before with 
aa arbitralily high weight value. If V is satisfactory, 
the design is complete. 

External drag is estimated by 

D = Cj) A 

where A is the external surface area of the nacelle. A is 
determined by dividing the nacelle in to three sections, 
the forebody from the nose to the point of tangency with 
a line parallel to the centerline at maximum diameter, a 
cylinderical midbody, and a tail section which is assumed 
to be symmetrical with the forebody. The forebody and 
tail can be approximated as paraboloids of revolution 
and their combined area is 

Ap = (47 t/3)(D^/2)2 [m - 1/(1 - M)] 

M = [ 1 - 4(2L/D^)2] i (ref. 4). 

For the cylinderical midbody, 

= n(D^2)2 (Ln - 2L) 
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and then A = A + A 

P ^ 

The drag coefficient is from Hoerner (ref. 15) and 
is given by 



Cj, = Cj [ 1 t 1.5 {D„/Ljj)3/2 * 7(D„Am)3] . 

The Schoenherr friction coefficient, C^, is calculated 
by a function subroutine, CFS(R), in which the argument 
is Reynolds number. 

The weight estimate depends on materials and details 
of the structure which were not the concern of this pro- 
ject. Therefore, only a rough approximation is possible 
and the assumptions are certainly open to question. 

Internal ducting weight was specified to be propor- 
tional to the cross sectional area and the length such 
that W = 15.07 (Area) (Length) in pounds 
where area is in square feet and length in feet (ref. 1). 

The weight of the outer envelope is estimated by 
assuming that it is proportional to the surface area, 
the maximum wall thickness, and the specific weight. 

A maximum wall thickness of 0.05Djjj was used to be con- 
sistant with the diffuser exit diameter assumption. 

The outer envelope of the nacelle, of course, is not 
solid but is made of inner and outer skins plus stiff- 
ening structure. The specific weight depends on the 
material which has not been specified. Two extremes of 
the weight spectrum which might be used for such a 
structure are aluminum and steel. Some material with 
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a weight between these two is more likely and a material 
density of one half that of steel was used in the sub- 
routine. The actual weight in water may very well be 
close to zero if the voids are filled with oil or pitch. 

The proportionality constant, a measure of the vol- 
ume of the material in the wall compared to the volume 
of the wall itself, a packing factor, will be a small 
decimal on the order of 0.01 to 0.1. 0.03 was used in 

the subroutine. 

The resulting weight was increased by ten percent 
to allow for the auxiliary inlet mechanism. 

At this point the subroutine’s job is completed and 
it returns to the calling program with the following 
quantities: 

a. the head loss summed from the free stream to the 
diffuser exit for both speeds, in feet, 

b. the drag of two nacelles for both speeds, in pounds, 

c. the exit area of two diffusers (one for each 
nacelle), in square feet, 

d. the weight of two nacelles, in pounds, 

e. the location of the vertical center of gravity of 
the nacelles below the keel of the craft, in feet, 
and 

f. the location of the longitudinal center of gravity 
of the nacelles forward of the transom, in feet. 
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Chapter 7« 

OFF-DESIGN PERFORMANCE 

Prediction of off-design performance may be carried 
out by SUBROUTINE NACEL only after a design has been com- 
pleted and returned to the calling program. The calling 
program must reset the value of a control variable named 
ISTRT from 1 which indicates design, to 3 indicating 

evaluation, and the value of a control variable named 

-9 

DELTA which has a finite value greater than 1 x 10 
during design to 0. which prevents rejection and return 
when cavitation limits are violated. 

The same initial values are known as when the design 
process started, and in addition, all values established 
during the design cycle are now stored for use during 
evaluation. 

The free stream conditions, total pressure, dynamic 
pressure and cavitation number, are calcuated for each 
speed to be evaluated. 

The subroutine now proceeds through the same set of 
calculations as for design except that the auxiliary in- 
let, and diffuser sizing sections are omitted. Cavitation 
checks are made in the same manner, but the existing 
parameters are not allowed to be rejected. Instead of 
rejection and return a variable named CAV(I,J) has its 
value changed from 0. to 1. to indicate that cavitation 

is occuring at the evaluation speed indicated by I and at 
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the location indicated by J. J=1 on the outside of the 
forebody, 2 inside the lip, and 3 in the elbow. 

On completion of an evaluation cycle the subroutine 
returns the value of the head loss in feet and the drag in 
pounds for two nacelles for each of the off-design speeds. 
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Chapter S. 

CONCLUSIONS AND RECOMMENDATIONS 

Results of the subroutine’s use with the optimiz- 
ation program are discussed in reference 1. In addition, 
to test it separately, the subroutine was called by a 
test program which stepped through values of from 

0.4 to 0.9 and values of cruise ^±/^q from 0.5 to 1.0 at 
three cruise speeds, 40, 45 and 50 knots. Take-off speed 
was assumed to be one half cruise speed. The required 
flow rate at cruise was arbitrarily set at 600 cubic 
feet per second with the flow rate at take-off ten per 
cent greater. The angle of attack was set at two degrees 
and depth of submergence at five feet. 

Ram pressure recovery at the diffuser exit and cav- 
itation limits as calculated from the above input are 
displayed in Figures IS through 21. The ram pressure 
recovery curves are unique to the particular flow rate 
combination assumed for this calculation and are present- 
ed as samples. The cavitation limits would be the same 
for any flow rate at the same angle of attack and depth 
of submergence. 

Several explicit and implicit assumptions were made 
in this study to simplify it or to try to bound a real- 
istic value. Further investigation might prove useful 
in any or all of the following areas: 

1. Very little information on bell diffusers that 
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could be directly applied in this study was found. 

Conical diffuser performance was assumed to be a lower 
bound on that expected from a bell. Information on bell 
diffusers similar to that in reference 11 on other shapes 
would allow realistic prediction of their performance. 

2. The pressure coefficient information used from 
reference 3 was calculated rather than measured. Although 
the predicted values for simple cases (i.e., a sphere) 
was compared with experimental results, no mention is 
made of comparison of calculated and experimental data 
for. a more complex shape such as an inlet nacelle. 

3. The drag calculation was made with the implicit 

assumption that the drag of a nacelle with a nose inlet 
and flow into that inlet is very close to the drag of 
that same body with no nose opening or inflow. Refer- 
ences 15, l6 and 17 support this assumption, but only a 
small number of cases with low are examined. There 

is some indication in reference l6 that the range of 
Vj^/V^ for which this assumption holds decreases as 
increases. 

4. The effect of proximity to the free surface on 
drag was not included. For overall system performance 
comparison, the omission of this effect would probably 
not invalidate the comparison but there could be a slight 
bias in favor of the smaller nacelle in the case where 

there is a significant difference in size. 
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In its present form, SUBROUTINE NACELLE requires 
interpolation in data tables several times. Some of these 
tables could be approximated by functions to speed up the 
operation and reduce computer costs. The data for which 
this could be carried out include Table 12 and Figure 17. 
The data in Tables 5 through 16 presents a more difficult 
problem because the relationship is not single valued in 
all cases. Perhaps separating it into two or more regions 
with a different function applicable in each region offers 
a solution. Some unsuccessful attempts were made to ex- 
amine possibilities for the large data bank. 

One other simplification which would decrease com- 
puter time would be accomplished by elimination of the 
calculation of diffuser length. One of the results of 
the data run mentioned above was that the diffuser length 
calculated by iteration in the subroutine was nearly 
always approximately ninety per cent of the mean of the 
maximum and minimum allowed values. The equation 

I^d “ 0»9(L(jniin "*■ ^dmax^/^ 

could be used to completely eliminate the iteration. 
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Table 1. Geometric Configuration of Kuchemann and Webber 
Circular Intakes Designated Class A, B and C> 

^ Dimensions of all three classes are defined by the 



following equations (adapted from 


I'Gf • 7 , p« 79)» 


Do/Di = 


Kl 




Dm/2L = 


K2 - K 3 (Ai/Am)3 




Di/2Li 


= I/K 4 (Ki - 1). 




Table l.a Values of the Constants. 




Class 


Ki K2 


K 3 


A 


1.15 0.2 


12.5 1.5 


B 


1.10 0.2 


12.5 1.3 


C 


1.03 0.2 


7.5 1.0 






(ref 


Table l.b Values of the Primary Ratios. 


Di/Dm 


L/D 


m 




Class A and B 


Class C 


0.90 


0.0732 


0.1195 


O.BO 


0.1439 


0.231 


0.70 


0.309 


0.452 


0.60 


0.639 


0.910 


0.50 


1.263 


1.575 


0.40 


2.00 


2.175 


0.20 


2.490 


2.494 



*(See Fig, 3 for illustration of dimension symbols.) 
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Table 2, Geometric Configuration of the Lockheed Intakes 



Nacelle 


*Do/Dm 


L/Li 


Di/Dm 


L/Dm 


A-40-550 


0.40 


9.0 


0.333 


2.20 


A-80-40 


o.go 


9.0 


0.77 


0.32 


A-50-2g0 


0.50 


9.0 


0.426 


1.40 


A-75-60 


0.75 


9.0 


0.713 


0.45 


A-70-90 


0.70 


9.0 


0.656 


0.63 


A-55-230 


0.55 


9.0 


O.4S3 


1.54 


A-60-135 


0.60 


9.0 


0.541 


O.gl 


A-45-400 


0.45 


9.0 


0.369 


1.60 


A-60-200 


0.60 


9.0 


0.541 


1.20 



(ref. 3) 

*(See Fig. 3 for illustration of dimension symbols.) 
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Table 3* Variation of Minimum Pressure Coefficient with 



Velocity Ratio 
A -40-550 
Internal Flow: 


and Angle 


of Attack for Lockheed 
Minimum Cp 


Nacelle 


Vi/Vo 


a=0 


20 


4® 


6° 


0.70 


0.420 


0.415 


0.400 


0.390 


O.SO 


0.260 


0.250 


0.240 


0.220 


0.90 


0.060 


0.050 


0.045 


0.030 


1.05 


-0.260 


-0.255 


-0.265 


-0.550 


1.15 


-0.500 


-0.510 


-0.550 


-1.040 


1.25 

External Flow: 


-0.760 


-0.765 


-1.050 


-1.630 


0.70 


-0.090 


-0.110 


-0.250 


-0.440 


0.60 


- 0.090 


-0.090 


-0.160 


- 0.320 


0.90 


- 0.090 


-0.090 


-0.120 


-0.230 


1.05 


-0.090 


-0.090 


-0.115 


-0.155 


1.15 


-0.090 


-0.090 


-0.115 


-0.145 


1.25 


- 0.090 


-0.090 


-0.115 


-0.140 



(ref. 3) 
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Table 4 * Variation of Minimum Pressure Coefficient with 
Velocity Ratio and Angle of Attack for Lockheed Nacelle 
A-gO -40 
Internal Flow: 



Minimum Cp 



Vi/Vo 


cxro 


2 ° 


4 ° 


6 ° 


0.70 


0.200 


0.050 






O.SO 


- 0.100 


- 0.235 






0.90 


-0.465 


-0.700 






1.05 


-1.160 


-1.425 






1.15 


- 1 . 7^0 


- 2.220 






1.25 


-2.500 








External Flow; 










0.70 


- 0.435 


-0.4^0 


-0 . 540 


- 0.60 


o.go 


-0.420 


-0.460 


-0.500 


-0.560 


0.90 


-0.405 


-0.450 


-0.490 


-0.515 


1.05 


-0.375 


-0.415 


-0.440 


- 0.475 


1.15 


-0.365 


-0.400 


-0.425 


-0.460 


1.25 


-0.350 


-0.375 


-0.410 


—0 • 440 



(ref .3) 
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Table 5» Variation of Minimum Pressure Coefficient with 
Velocity Ratio and Angle of Attack for Lockheed Nacelle 



Minimum 



A-50-2CO. 

Internal Flow: 

Vi/Vo 

0.70 

0.^0 

0.90 

1.05 

1.15 

1.25 

External Flow; 
0.70 
0.80 
0.90 
1.05 
1.15 
1.25 



ot=0 


20 


0.380 


0.355 


0.200 


0.160 


-0.015 


-0.050 


-0.370 


-0.420 


-0.640 


-0.700 


-0.940 


-1.00 


-0.150 


-0.160 


-0.150 


-0.155 


-0.150 


-0.155 


-0.150 


-0.155 


-0.150 


-0.155 


-0.150 


-0.155 



Cp 

4 ° 6 ^ 
0.320 
0.120 
- 0.100 
-0.485 
-0.770 
- 1.080 

-0.240 

- 0.200 

-0.190 

-0.175 

-0.175 

-0.175 

(ref. 3) 
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Table 6« Variation of Minimum Pressure Coefficient with 
Velocity Ratio and Angle of Attack for Lockheed Nacelle 
A-75-60 
Internal Flow: 

Minimum Cp 



Vi/Vo 

0.70 

O.gO 

0.90 

1.05 

1.15 

1.25 

External Flow: 
0.70 
0.^0 
0.90 
1.05 
1.15 
1.25 



oc =0 


2° 


0.250 


0.130 


-0.020 


-0.175 


-0.340 


-0.530 


-0.930 


-1.230 


01.440 


-I.S15 


-2.050 


-^-2.50 


-0.360 


-0.390 


-0.350 


-0.375 


-0.340 


-0.360 


-0.320 


-0.345 


-0.310 


-0.335 


-0.300 


-0.325 



40 6° 

- 0.020 
-0.350 
-0.7^0 
- 1.600 
-2.370 
^-2.50 

-0.445 

-0.410 

-0.395 

-0.370 

-0.355 

-0.345 



(ref. 3) 
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Table ?• Variation of Minimum Pressure Coefficient with 



velocity Ratio 
A-70-90. 
Internal Flow: 


and Angle 


of Attack for Lockheed 
Minimum Cp 


Nacelle 


Vi/Vo 


II 

o 


20 


40 


6° 


0.70 


0.305 


0.215 


0.100 




o.go 


0.070 


-0.040 


-0.1^0 




0.90 


-0.210 


-0.350 


-0.535 




1.05 


-0.710 


-0.930 


-1.225 




1.15 


-1.130 


-1.420 


-1.920 




1.25 

External Flow: 


-1.630 


-2.120 


^ -2.50 




0.70 


-0.2S0 


-0.310 


-0.375 


-0.455 


0.^0 


-0.270 


-0.295 


-0.335 


-0.405 


0.90 


-0.265 


-0.2^5 


-0.315 


-0.365 


1.05 


-0.255 


-0.270 


-0.300 


-0.330 


1.15 


-0.250 


-0.265 


-0.290 


-0.320 


1.25 


-0.240 


-0.260 


-0.2^0 


-0.310 



(ref. 3) 
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Table Variation of Minimum Pressure Coefficient with 
Velocity Ratio and Angle of Attack for Lockheed Nacelle 



Minimum 



A-55-2gO.. 
Internal Flow: 

Vi/Vo 

0.70 

O.SO 

0.90 

1.05 

1.15 

1.25 

External Flow: 
0.70 
O.SO 

0.90 

1.05 

1.15 

1.25 



<x =0 


2° 


0.415 


0.3S5 


0.235 


0.200 


0.030 


-0.010 


-0.330 


-0.375 


-0.595 


-0.650 


-o.aso 


-0.950 


-0.125 


-0.225 


-0.125 


-0.145 


-0.125 


-0.135 


-0.125 


-0.130 


-0.125 


-0.130 


-0.125 


-0.130 



Cp 



4® 


60 


0.350 




0.160 




-0.070 




-0.440 




-0.7B0 




-1.3^0 




-0.450 


-0.700 


-0.270 


-0.500 


-0.1^5 


-0.320 


-0.170 


-0.205 


-0.150 


-0.190 


-0.150 


-0.170 



(ref. 3) 
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Table 9» Variation of Minimum Pressure Coefficient v/ith 
Velocity Ratio and Angle of Attack for Lockheed Nacelle 
A-60-135 
Internal Flow: 

Minimum Cp 



Vi/Vo 

0.70 

o.ao 

0.90 

1.05 

1.15 

1.25 

External Flow: 
0.70 
0.80 
0.90 
1.05 
1.15 
1.25 



cc =0 


20 


0.300 


0.230 


0.080 


-0.010 


-0.175 


-0.280 


-0.620 


-0.755 


-0.965 


-1.120 


-1.340 


-1.535 


-0.255 


-0.280 


-0.255 


-0.275 


-0.250 


-0.270 


-0.245 


-0.265 


-0.240 


-0.265 


-0.240 


-0.260 



4° 


60 


0.140 


0.040 


-0.115 


-0.235 


-0.405 


-0.560 


-0.910 


-1.135 


-1.330 


-1.620 


-1.80 


-2.230 


-0.315 


-0.380 


-0.300 


-0.350 


-0.295 


-0.330 


-0.280 


-0.310 


-0.280 


-0.300 


-0.275 


-0.295 




(ref. 3) 
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Table 10. Variation of Minimum Pressure Coefficient with 
Velocity Ratio and Angle of Attack for Lockheed Nacelle 
A-45-400 
Internal Flov/; 

Minimum Cp 



Vi/Vo 

0.70 

0.^0 

0.90 

1.05 

1.15 

1.25 

External Flow; 
0.70 
O.SO 
0.90 
1.05 
1.15 
1.25 



ot=0 


20 


0.400 


0.3^5 


0.230 


0.205 


0.025 


0.000 


-0.320 


-0.340 


-0.570 


-0.600 


-0.^50 


-o.aso 


-0.120 


-0.140 


-0.120 


-0.130 


-0.120 


-0.130 


-0.120 


-0.130 


-0.120 


-0.130 


-0.120 


-0.130 



40 60 

0.370 

0.1^5 

- 0.020 

-0.365 

-0.630 

-0.915 

-0.310 

- 0.210 

“0.155 

-0.150 

-0.150 

-0.150 

(ref. 3) 
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Table 11. Variation of Minimum Pressure Coefficient with 
Velocity Ratio and Angle of Attack for Lockheed Nacelle 
A-60-200 
Internal Flow: 



Minimum Cp 



Vi/Vo 


ot=o 


2° 


4° 


60 


0.70 


0.400 


0.360 


0.300 




0.^0 


0.205 


0.150 


0.075 




0.90 


-0.015 


-O.OSO 


-0.165 




1.05 


-0.400 


-0.4^5 


-0.620 




1.15 


-0.6S5 


-0.^00 


-1.140 




1.25 


-1.000 


-1.265 


-1.^90 




External Flow: 










0.70 


-0.160 


-0.200 


-0.3^5 


-0.660 


O.SO 


-0.160 


-0.175 


-0.240 


-0.420 


0.90 


-0.160 


-0.160 


-0.210 


-0.2S0 


1.05 


-0.160 


-0.160 


-0.1^5 


-0.230 


1.15 


-0.160 


-0.160 


-O.lSO 


-0.215 


1.25 


-0.160 


-0.160 


-O.lSO 


-0.210 



(ref. 3) 
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Table 12, Pressure Loss Coefficient for an Elliptical 
Profile Lip 



Vi/Vo 


AP/qo 


0.6 


0.025 


O.S 


0.033 


1.0 


0.043 


1.2 


0.0S5 


1.4 


0.130 


1.6 


0.173 



Reynold's number^S X 10^ 

(Adapted from Fig, 74, Report No, 5, ref. 3) 
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Figure 3» Nacelle Geometry 



I 



I 



2.4 




Figure 4* L/Dm vs. Di/Dm for the Lockheed intakes and 
Kuchemann and Weber intakes Classes A, B and C. 
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Figure 5» Minimum Cp vs. L/Dm in the external flow at Vi/Vo 
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,5. 1.0 L/Dm 1.5 2.0 

Figure 6. Minimum Cp vs. L/Dm in the external flow at Vi/Vo = .50 
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Figure S, Minimum Cp vs. L/Dm in the external flow at Vi/Vo = 1.05 
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Figure 9» Minimum Cp vs. L/Dm in the external flow at Vi/Vo = 1.15 
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Figure 10. Minimum Cp vs. L/Dm in the external flow at Vi/Vo = 1.25- 
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Figure 11. Minimum Cp vs. L/Dm in the internal flow at Vi/Vi 
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Figure 12. Minimum Cp vs. L/Dm in the internal flow at 
Vi/Vo = .go . 
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Figure 13. Minimum Cp vs. L/Dm in the internal flow at 
Vi/Vo = .90 
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Vi/Vo =1.05 
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Figure 15. Minimum Cp vs. L/Dm in the internal flow at 
Vi/Vo =1.15 
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Figure 16. Minimum Cp vs. L/Dm in the internal flow at 
Vi/Vo = 1.25 
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Figure 17. Ram Pressure Recovery vs. for the elliptic- 

al profile lip (Blackaby and Watson IBE) . 
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Figure IS. Ram Pressure Recovery at the diffuser exit vs. 
Vj^/V^ for a 40 knot cruise speed 



- 73 - 




Figure 19 . Ram Pressure Recovery at the diffuser exit vs. 
V^/V^ for a 45 knot cruise speed 
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Figure 20. Ram Pressure Recovery at the diffuser exit vs. 
for a 50 knot cruise speed 
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Figure 21. Cavitation limits as a function of D./D vs. 
^ i' ra 

V^/V^ at cruise speeds of 40, 45 and 50 knots 



- 76 - 



Appendix A. CAVITATION CRITERA 

The total pressure head at any point in the system can 
be determined with respect to the free stream head by trac- 
ing along a stream line from out in the free stream to 
the point of interest. 



H - Hq - 

or P = Pq- Pioss terms of pressure. 

Using Bernoulli’s equation and allowing the static 
head term to reflect all changes in elevation or 
submergence, 



h„ - v2/2g = h - v2/2g - 

or . i/>v2 = p * 

The cavitation number, cr t can be defined as 

(Pref. - 

Where Pj.^^ and V^^^ are the reference pressure and 
velocity to which the cavitation number is related and 
p is the local static pressure at the point of interest. 

If the free stream conditions are used for reference, 
then Pq and Vq are the reference pressure and velocity. 

In this case 



P„ - P = i V2 - i V2 . 

so = (Po - puml = v2/v2 - 1 ^ PiossA/^v2 

Cavitation is assumed to occur when the local static 
pressure is equal to the vapor pressure of water. The 
incipient free stream cavitation number is defined by 

t^oi : (Po - Pyj/i/vZ P v2)-l * Pi„33/i/>v2. 
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This definition allows the prediction of cavitation using 
the free stream incipient cavitation number if the losses 
are known. 

Consider the situation in Figure (3)» The diffuser 
entry is designated station 1, and the diffuser exit 
station 2, to remain consistant with the rest of the text. 
The turning vane is not intended to be a particular type 
or design, but is expected to be well designed and sized 
for the elbow and diffuser with which it is used. The 
sketch is conceptual, only. The subscript tv will designat 
conditions on the turning vane. 

To determine the cavitation number on the turning vane 
conditions at station 2 can be used and, in turn, related 
to the free stream. As before, 

^tv “ ^2 ” ^loss 

2 — >tv 

If the turning vane starts near station 2, then any 
loss between there and the point of lowest pressure can 
be neglected as it will be very close to the leading edge. 

hv = ^2 = Ptv * = P2 

The cavitation number on the turning vane is, 

<5tv = <P2 - 

When p^^ = p^, cavitation starts, so the incipient 
cavitation number is 

^tvi " (P2 “ 

(T is a characteristic number determined from 
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experiments and is therefore known (ref. 28). The con- 
ditions at station 2 which will cause cavitation on the 
turning vane can be determined. Solving the equation 
above for p^, 

P2 = <^tvi ^’'2 * Pv> 

but P 2 can be determined in terms of the free stream 
conditions, which are known. 

Pa = Po ^ - ho33 



^vi ^^^2 + Pv = Po ■ ^'^^2 ■ 



loss 



^loss this expression is the sum of all losses from the 
free stream to station 2. ^2 only variable re- 

maining which indicates the conditions at station 2, so 
2 

solving for V^, 

Vo (1 + CT .) = f(p - P )/s/>V^ + 1 - P /s/’V^l 
2 "^tvi L'*"o o loss oJ o 

i 

The first term in the brackets is the free stream 
incipient cavitation number. The value of which gives 
rise to incipient cavitation on the turning vane is 

loss 



''2 = I 



<7-oi * 1 - Pi 



]i vy (1 4 



- 79 - 



Appendix B. GEOMETRIC DERIVATIONS 

There are three variable lengths in the nacelle which 
are not a function of the forebody geometry and which must 
be calculated in each design. They are the diffuser lengtn , 
L^, the additional length due to the auxiliary inlet, L^, 
and a length after the diffuser exit containing part of 
the elbow which directs the flow upward into the strut, the 
strut inlet length. 

1. The Strut Inlet Length 

Assume that the thickness to chord ratio of the strut 
all along its length from the nacelle to the free surface 
is determined by the criterion 

^mAo "" [(1 + l] A-15 (ref. 1,4) 

where cr± is the local free stream incipient cavitation 
number. Also assume an elliptical cross section for the 
strut and a rectangular cross section for the internal 
ducting (ref. 1). 

The respective areas are 

^strut “ (TT^/4) tjjj Cq 

^duct ~ i ^m *^o (ref. 1) . 

Aduct '^^st be the same as the exit area of the diffuser, 

A2. 

To allow for structure and fairing, assume that the 
desired additional nacelle length is close to the value of 
the strut chord, Cq. 

For CT - 0.3, tj^ = 0.12 c^, 2 A 2 = 0.12 so that 



- 80 - 



Co = 3.545 Dj . 

2. Length Due to the Auxiliary Inlet 

If the auxiliary inlet has the configuration shovm 
in Figure 3, the inlet area is a section of a right circ- 
ular cone with a minimum diameter and slant height X. 

The angle of inclination from the centerline of that 
cone section, shown in the figure as is assumed to be 
formed between the centerline and a line passing through the 
nose leading edge to the point of tangency of the forebody 
at maximum diameter, 

Sx ' tan-1 [(D_^ - Di)/2L] 

The mean diameter of the conical section is then 
D = Dj_ + X sin ©X 

and its area is 

Acone = D X . 

The area of the inlet will actually be some percentage 
of the area above because some allowance must be made for 
structure. 

hux - K nD X 
K = 0.© is used in the subroutine. 

Finally, solving for X in the above expression, 

X = (0.5/sin 0^) [(Di^ + i.273 sin ©^/K)^- Di] , 

and = X/cos ©^ . 

3. Diffuser Length 

The length of the diffuser, L^, becomes a factor in 
the determination of nacelle length when the sum of the 
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internal lengths exceeds the minimum length, The 

excess length is chargeable to the diffuser because all the 
other internal length requirements are fixed by other con- 
siderations. The difference in length is 

= ^Lip + Lx ^ 3.545 D2 - 5.5 Dtn . 

The change in drag is given by 
A Drag = 

because Cq is essentially a constant for small changes in 

The change in length can be considered to be all in a 
cylinderical section of diameter 

A. S - 7T Djjj A 

The increase in power due to the increase in drag is 
A Powerj^ = Cp A Lfj . 

The power to overcome the diffuser losses is found by 
APower2 = (CK^ + . 

The total increase in power required to operate the 
diffuser is then 

A Power = APower]^ + A Power2 . 

In the subroutine, the minimum value of A Power is 
determined by iteration on 
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Appendix C. SUBROUTINE LISTING 



SUBROUTINE NAGEL and its flow diagram listed in the 
following pages are in the form for v;hich the rationale is 
described in chapters 6 and 7. Similar listings for 
FUNCTION TABLE, which was not programmed by the author but 
was furnished for common ues with reference 1, follow 
NACEL« FUNCTION CFS is also shown for completeness but no 
flow diagram was generated for it because of its simplicity. 

The FORTRAN names below, used for various variables and 
constants in the subroutine, are described as an aid in 
understanding the logic and specific calculations. The 
symbols are listed in approximate order of first occurance 
beginning with the first executable statement. 

ZK - A constant decimal less than one which in- 



HA 



SPO 



HS 



dicates how much of the annulus occupied by 
the auxiliary inlet is available for opening. 
Static pressure in the free stream in p.s.f. 
Depth of submergence in feet. 

Pressure head of the atmosphere in feet of 
water. 



RHOW - 



Density of water in slugs per cubic foot. 
Acceleration of gravity, 32.174 feet/sec? 



G 



PVP - Vapor pressure of water in p.s.f. 

PV ~ Vapor pressure of water in feet of water. 

VELR{I) - Inlet velocity ratio at speed I. 1=1 



indicates cruise. 1=2 indicates take-off 
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VI(I) - 
VO(I) - 
SIGTV - 



JNUMB - 
ISTRT - 



NUMB - 

QO(I) - 

SIGI(I) - 
PTO(I) - 
TRIM(I) - 
CPEX - 
CDUMX - 
CDUMY - 
CPEXT(J,K 



I greater than 2 indicates a performance 
evaluation speed. 

Inlet velocity in feet/sec. 

Free stream velocity in feet/sec. 

Incipient cavitation number on the turning 
vanes referenced to diffuser exit pressure 
and velocity. 

A counter. 

A control variable indicating whether the 
current calculation is part of the design 
process (ISTRT =1), or part of the eval- 
uation process (ISTRT = 3)» 

A control variable indicating the total 
number of speeds to be considered. 

Free stream dynamic pressure at speed I in 
p.s.f . 

Free stream incipient cavitation number. 

Free stream stagnation pressure in p.s.f. 

Angle of attack in degrees. 

External minimum pressure coefficient. 

A dummy array generated during interpolation 

A dummy array generated during interpolation 

,1) - Tabulated minimum pressure coefficients 

in the external flow at the Jth velocity 

ratio, Kth length to diameter ratio, and Ith 

angle of attack. 
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CPINT(J,K,I) - Tabulated minimum pressure coefficients 
in the internal flow, 

VRT - Tabulated array of velocity ratios. 

IL, KL, liL, NL - Control variables indicating the 
order of interpolation desired. 

XDT - Tabulated array of angle of attack. 

XDTT(I) - A dummy array generated during interpolation. 

DIDMT - Tabulated array of Di/Dj^. 

DIDMX - The maximum value of D^/Djj, for a V^/Vq as 
determined from interpolation in the data 
table. 

XD - A working variable for length to diameter 

ratio. 

WGTS(1,1) - Wet weight of the nacelle in pounds. 

DIDM - Current value of 

QI - Flow rate at cruise through one inlet in 

c . f . s . 

AI - Nose inlet area in feet squared. 

DI - Current value of in feet. 

DM - Current value of Djjj in feet, 

ELEXT - Length of the forebody in feet. 

VRTT(I) - A dummy velocity ratio array generated during 
interpolation. 

VRTEX(I) - A dummy velocity ratio array generated during 
interpolation . 

VRMX(I) - Maximum allowed velocity ratio for a given 
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QIN - 

QC - 
QAUX - 

KDEX - 
VI2 - 
PRLT - 

PRLT2 - 

PTI - 

SPI - 
AIAUX - 
VIAUX - 
PTAUX - 

PRAUX - 

DYP - 

PC - 



D^/Djjj as determined from the tables. 

Working variable for flow rate through the 
nose inlet at take-off in c.f.s. 

Required flow rate at take-off in c.f.s. 

Flow rate through the auxiliary inlet in 
c . f . s . 

A counter. 

Working variable for VI(2). 

Tabulated array of pressure recovery coeffic- 
ients for the lip. 

Lip pressure recovery coefficient at take- 
off. 

Stagnation pressure in the internal flow in 

p.s.f . 

Static pressure in the internal flow in p.s.f. 
Area of the auxiliary inlet in feet squared. 
Auxiliary inlet velocity ratio. 

Stagnation pressure inside the auxiliary 
inlet in p.s.f. 

Pressure recovery coefficient of the aux- 
iliary inlet. 

Dynamic pressure inside the auxiliary inlet 
in p.s.f. 

Stagnation pressure in the combined flows 

from the auxiliary and nose inlets in p.s.f. 

Inlet losses at speed I in p.s.f. 
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DLIP(I) 



QDIF - 


Flow rate through the diffuser in c.f.s. 


ELENT - 


Interior length of the lip in feet. 


PHI - 


The angle labeled 6^ in Figure 3» 


ELAUX - 


Interior length due to the auxiliary inlet 
in feet. 


ELMAX - 


Maximiim permissible length of the diffuser 
in feet. 


ELMIN - 


Minimum permissible length of the diffuser 
in feet. 


QDIF - 


Flow rate through the diffuser in c.f.s. 


EL - 


Current value of the diffuser length in feet 


DEL - 


Change in diffuser length in feet. 


ELN - 


Length of the nacelle in feet. 


ELD, ELL 


- Working variables for ELN. 


DDM - 


Mean diameter of the diffuser in feet. 


XKT - 


Sudden expansion loss coefficient, K^. 


REL - 


Reynolds number based on length. 


RED - 


Reynolds number based on diameter. 


DL - 


Ratio, Din/Lj^. 


CDRG - 


Drag coefficient. 


ANGL - 


Deffuser half angle of expansion, labeled 
6 in Figure 3 • 


POW - 


Current value of the power charged to the 
diffuser in ft. lbs. /sec. 


POWI - 


Previous value of POW. 


ELDIF- 


Final value of diffuser length in feet. 
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REND - 
RENL - 
DDIF - 
FLOSS - 
VAOUT - 
VCRIT - 



VI-lAX - 
DELH(I,1 

PRC(I) - 

CD(I) - 
AEXN - 
CPOD(I) 
AREA(l) 

CGS(1,1) 

CGS(2,1) 

CGS(3,1) 

CGS(4,1) 



Same as RED. 

Same as REL. 

Diffuser loss in p.s.f. 

Sum of diffuser, pipe and lip losses in p.s.f. 
Average diffuser exit velocity in ft ./sec. 
Critical diffuser exit velocity at which 
cavitation occurs on the turning vanes in 
ft. /sec. 

Maximum diffuser exit velocity in ft. /sec. 

) - PLOSS converted to head in feet of water 
at speed I. 

Ram pressure recovery coefficient at the dif- 
fuser exit. 

Drag coefficient. 

Wetted area of the nacelle in feet squared. 

- Drag of two nacelles in pounds. 

- Diffuser exit area of two diffusers in feet 
squared. 

- Distance of the vertical center of gravity 
of the nacelles from the keel in feet. 

- Distance of the longitudinal center of 
gravity of the nacelles from the transom in 
feet. 

- Distance of the VCG of the contained water 
from the keel in feet. 

- Distance of the LCG of the contained water 
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RHOD - 



from the transom in feet. 

Density of steel in slugs per cubic foot. 
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VRMAX (1 ) = TAaLEI ALFAT, VRTFX ,TR I M( 1 ) , ML ) 
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Appendix D. ALTERNATE DESIGN RATIONALE USING ONLY VELOCITY 
RATIO AS A SYSTEM PARAMETER 



Another .approach to the design of the inlet nacelle 
can be made using the maximum allowable at a given 

Vj^/V^ in stead of allowing it to vary as a parameter. The 
design proceeds as previously described except that the 
diameter ratio determined from interpolation in the data 
table is used to calculate the nacelle geometry. If, after 
the diffuser is sized, the exit velocity is too high, then 
instead of rejecting that set of parameters, a new value 
of Djjj based on the maximum allowed exit velocity is cal- 
culated and the subroutine goes back to the beginning of 
the design phase where a new design begins with this value* 
Performance evaluation proceeds exactly as before. 

The following pages contain a listing of the subroutine 
in this form with a flow diagram. 
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DATA ALFAT/0. 0,2. 0,4. 0,6.0/, 

XDT/0. 25, 0.50, 0.75,1 .0, 1, 2 5,1.5 0, 1.75 , 2. 0 , 2. 2 5 ,2 . 5 / , D I DMT /O. 8 10 
0.695, 0.627,0.577,0.533, 0 . -^90 , 0. 4 39 , 0 . 365 , 0. 260, 0. 090/ 

, VRT/. 7, .8, ,9, 1.05, 1.15, 1.25/ 



DATA AA/ 
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Appendix E. ASSUI-ffiD CONSTANTS AND LIMITS 

It was necessary to make numerous assumptions and 
impose limits on various quantities in the subroutine. Some 
of the assumptions were somewhat arbitrary and may be mod- 
ified whenever good reasons to do so are discovered. 

Some of those assumptions and limits are discussed 
here in the approximate order in which they are encountered 
in SUBROUTINE NACEL. 

The first assumption that should be mentioned 
is that everything in the model is based on a nacelle with 
a circular cross section. No attempt was made to account 
for other shapes. 

The decimal constant, ZK, which determines the 
amount of available inlet area in the auxiliary inlet ring 
has an arbitrary value of O.S. The real case may vary 
widely from this depending on the details of the design. 

The incipient cavitation number on the turning 
vanes, SIGTV, has a value of 0.4. This value was accepted 
as a reasonable one for the speed range of interest. No 
examples with which to compare were found in the literature. 

The angle of attack at cruise, TRIM(l), has been 
artificially limited to a maximum of three degrees because 
the external pressure coefficients are not single valued 
with respect to L/Djjj or Di/Dj^ for ck ^ 4 °. The highest 
tabulated value for angle of attack is 6° and extrapolation 
would probably not be valid. 
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The range of tabulated velocity ratios is from 
0.7 to 1 . 25 . Extrapolation is allowed in the subroutine 
but the results could contain large errors if extrapolated 
very far. 

The largest value of L/Dj^ in the tables is 2.50. 
This is also allowed to be extrapolated and the same remarks 
apply as for the velocity ratios. 

is limited to O. 9 O, which assumes a wall 
thickness allowance of 0.05Djjj. This value is arbitrary and 
may be modified as desired. 

The pressure loss factor for the auxiliary inlet 
is assumed to be 0.2 which contains an underlying assumption 
that its inlet velocity ratio is always near 1.0. 

Diffuser lengths are limited to those correspond- 
ing to double angles of 6° to 20®. See chapter 4 for an 
explanation of the limits. 

In the wetted surface calculation, the area of 
fairings into the foils and strut, as well as the missing 
area in way of the foils and strut are ignored. The nacelle 
is assumed to be a solid body of revolution. 

Interference drag between the nacelle and the 
other components is assumed to remain constant even though 
the nacelle size changes. No attempt is made to account 
for the effect of or the depth of submergence on drag. 

Chapter S contains recommendations concerning this subject. 
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The weight estimate, as discussed in chapter 6, 
is only a rough approximation with a series of assumptions 
as to the construction and m.aterial. 
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